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SUMMARY 

1 

Lifting-surface  calculations  for  a  particular  wing  at  Mach  number  0.8 
are  used  to  examine  the  aerodynamic  consequences  when  one  side  of  an  air¬ 
craft  enters  a  vertical  gust  before  the  other.  The  growth  of  lift,  pitch¬ 
ing  moment,  aileron  hinge  moment  and  the  behaviour  of  the  transient  rolling 
moment  are  analysed.  Their  dependence  on  the  inclination  of  the  gust  front 
and  the  length  of  a  (1-cosine)  gust  is  discussed. 

Some  idealized  calculations  of  control-surface  motion  to  suppress  the 
wing  root  bending  moment  or  to  achieve  undisturbed  flight  show  only  minor 
effects  of  gust  inclination  on  control  demand.  But  the  transient  rolling 
moment  from  asymmetric  entry  into  a  (1-cosine)  gust  within  the  spectrum  of 
current  airworthiness  requirements  is  large  enough  to  constitute  a  consider¬ 
able  handling  problem.  Ride  control  studies  to  eliminate  heave,  pitch  and 
roll  by  means  of  two  ailerons  and  an  auxiliary  tailplane  or  canard  reveal 
distinct  advantages  of  the  canard  as  an  active  control. 
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1  INTRODUCTION 

Airworthiness  requirements  for  an  aircraft  structure  include  the  ability  to  survive 
gust  loads  under  prescribed  conditions  with  an  ultimate  safety  factor  of  1.5.  Aircraft 
are  designed  for  discrete  gust  cases  involving  ramp  or  (1-cosine)  gusts  with  velocities 
approaching  70  ft/s  (21  m/s).  Horizontal  gusts  are  considered  in  order  to  ensure  adequate 
stressing  of  the  fin,  while  the  horizontal  lifting  surfaces  must  safely  bear  the  additional 
structural  loads  in  specified  vertical  gusts  in  the  absence  of  trimming  action  by  the 
pilot . 

It  is  sufficient  to  demonstrate  this  airworthiness  quality  for  the  aircraft  in  level 
flight  at  constant  speed  in  the  direction  normal  to  the  gust  front.  This  last  simplifica¬ 
tion  conveniently  maintains  spanwise  symmetry  in  aerodynamic  loading,  but  in  reality  one 
wing  will  enter  a  gust  before  the  other  and  the  aircraft  will  experience  rolling  moment 
as  well  as  pitching  moment  and  a  change  in  lift.  The  primary  purpose  of  the  present 
investigation  is  to  calculate  these  time-dependent  forces  for  different  flight  directions 
relative  to  the  gust  front  for  a  step  gust  and  also  for  (1 -cosine)  gusts. 

The  wing  to  be  considered  is  that  used  in  the  theoretical  study  of  active  control 
aerodynamics  in  Ref  1,  and  the  lifting-surface  calculations  will  again  be  limited  to  the 
particular  Mach  number  0.8.  The  planform  of  aspect  ratio  8  is  shown  in  Fig  I  and  includes 
port  and  starboard  ailerons,  which  will  be  used  as  independent  active  controls  to 
neutralize  the  rolling  moment  and  root  wing  bending  moment.  These  purely  illustrative 
results  will  give  some  indication  of  the  influence  of  asymmetry  of  gust  entry  on  control 
demand  in  structural  load  alleviation.  Pitch  control  will  be  achieved  by  means  of  an 
additional  tailplane  or  canard  surface,  so  that  the  three  degrees  of  freedom  in  heave, 
pitch  and  roll  can  be  considered  simultaneously.  The  achievement  of  perfect  ride  control 
through  a  gust  with  three  active  controls  is  an  entirely  aerodynamic  and  somewhat  academic 
exercise.  Nevertheless,  the  calculated  wing  root  bending  moment  will  serve  tp  illustrate 
how  a  typical  structural  characteristic  may  depend  on  asymmetric  gust  entry. 

2  GROWTH  OF  AERODYNAMIC  FORCES 

The  fundamentals  of  indicial  aerodynamics  in  linearized  compressible  flow  have  been 

2  3 

reviewed  by  Lomax  ,  whose  earliest  references  are  to  the  classical  work  of  Wagner  and 

..  U  ...  .  . 

Kussner  for  two-dimensional  incompressible  flow.  The  Wagner  function  relates  to  a 

single  abrupt  change  in  upwash  at  a  lifting  surface,  and  its  modem  equivalent  for  the 

various  aerodynamic  forces  due  to  a  step  change  in  aileron  angle  is  used  in  section  3.1. 

We  first  consider  the  counterpart  of  the  Kussner  function  for  the  forces  on  a  swept  wing 

as  it  enters  a  stationary  vertical  gust  of  constant  velocity  in  level  flight  at  an 

arbitrary  inclination  p  to  the  normal  to  the  gust  front. 

As  in  Ref  5,  approximate  estimates  of  lift,  pitching  moment  and  rolling  moment  due 

*  £ 

to  asymmetric  entry  into  a  step  gust  can  be  made  by  means  of  piston  theory  .  But  very 
few,  if  any,  such  calculations  appear  to  have  been  made  on  the  basis  of  lifting-surface 
theory.  Our  primary  objective  is  to  study  the  influence  of  the  flight  direction  p  on 
various  Kussner  functions  with  the  aid  of  the  subsonic  lifting-surface  theory  of  Davies^. 
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2 . 1  Oscillatory  calculations  with  spanwise  asymmetry 

Essential  asymmetry  in  lifting-surface  calculations  relatively  uncommon.  It  has 
arisen  in  the  study  of  acoustic  propagation  in  Ref  8,  which  is  focussed  on  the  spanwise 
loading  due  to  a  single  oscillating  aileron.  Antisymmetric  perturbations  of  a  symmetric 
mean  flow  are  normally  such  that  the  two  parts  are  independent,  but  there  is  a  case  of 
interaction  in  Ref  9  where  a  second-order  quantity,  the  yawing  moment  due  to  roll  rate, 
involves  products  of  symmetric  and  antisymmetric  parts  of  the  spanwise  loading. 

To  treat  asymmetric  gusi  entry,  we  first  consider  a  sinusoidal  gust  with  upwash 
velocity 

where  U  is  the  flight  speed  in  the  negative  x-direction,  y  is  to  starboa-d,  and 
Xq(iJj)  is  the  value  of  x  where  the  gust  front  intersects  the  centre-line  y  =  0  at 
the  moment  of  entry.  The  complex  upwash  angle  with  frequency  parameter  v  =  mc/U  is 


w(*,y)  -  co- 

^x  —  y  tan  t ]/ -  x Ai>) 

.1 

-  i  sin 

p 

x  -  y  tan  \(i  -  xQ(t) 

1 

u  cos 

L  5  J 

L  *  J 

which  is  split  into  its  symmetric  and  antisymmetric  parts 


(2-2) 


w(x,y)  -t-  w(x,-y) 
2U 


and 


w(x,y)  -  w(x,-y) 
2U 


(2-3) 


in  preparation  for  the  lifting-surface  calculations. 

Section  2.1  of  Ref  1  includes  some  comments  on  calculations  with  the  method  of 
Davies  .  Each  force  mode  is  written  as 


and  each  upwash  mode  as 


z  =  ftjz.  (x,y)eia)tJ 


w 

U 


ft 


w. (x,y) 
U 


(2-4) 


(2-5) 


All  these  quantities  are  usually  symmetric  or  antisymmetric  in  the  spanwise  direction. 
The  corresponding  lift  per  unit  area 


ft^pU2? .  (x,y)eiu)tj 


(2-6) 


with  stream  density  p  is  obtained  from  numerical  solution  of  an  integral  equation 
relating  w^  and  ,  whence  the  complex  generalized  forces  are  given  by  integrals 


a 


) 
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z^x.y)  2 ^ (x,y)dxdy 


(2-7) 


over  the  planform.  For  our  purposes  it  is  sufficient  to  take  integers  i  and  j  from 
1  to  5.  Modes  1  to  3  correspond  to 


z(  =  c  (lift) 

z_  =  x  -  x  (nose-down  pitching  moment) 
Z  3 


(2-8) 


z3  =  x  -  xh 


on  the  aileron  surfaces  and  zero  over  the  rest 
of  the  planform  (hinge  moment) 


and  the  corresponding  upwash  modes 


w.  dz.  iuz . 

— 1  =  — =L  +  _ 1 

U  3x  U 


(j  =  1,2,3) 


(2-9) 


represent  heave,  pitch  and  symmetric  aileron  motion  respectively.  It  is  convenient  to  fix 
the  origin  of  coordinates  where  the  pitching  axis  intersects  the  centre-line  (Fig  I),  so 
that  x^  =  0  .  The  subscript  4  refers  to  antisymmetric  conditions,  but  the  force  and 
upwash  modes  are  unrelated.  We  take 


z^  =  y  (rolling  moment) 

=  fy/lyl)w3  (antisymmetric  aileron  motion) 


(2-10) 


The  only  other  force  mode  to  be  used  is  the  root  bending  moment  given  by 


z,.  =  J y  (  (bending  moment) 


(2-1  1) 


and  the  integer  j  =  5  denotes  the  sinusoidal  gust,  so  that  w,.(x,y)  is  unconventionally 
asymmetric  in  y  unless  <p  =  0  and  is  given  by  equation  (2-2).  Thus  the  first  of 
expressions  (2-3)  is  used  for  w^  when  the  force  mode  is  symmetric  (i  =  1,2, 3, 5)  and  the 
second  one  is  used  for  the  rolling  moment  (i  =4).  To  obtain  the  separate  root  bending 
moments  on  the  port  and  starboard  wings,  we  suppose  that  if;  is  positive  and  the  starboard 
wing  with  y  positive  is  the  leading  side  relative  to  the  gust.  Let  B„  and  denote 

the  bending  moments  on  the  leading  and  trailing  (port)  sides  respectively.  The  rolling 
and  total  root  bending  moments  are  thus 


and 


whence 


B  "  B*.  +  Bt  - 

B(  =  icU2c3(Q55  ♦  Q45) 

Bt  =  JouV(Q55  -  Q45) 


(2-12) 

(2-13) 


(2-14) 
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At  the  outset  of  the  lifting-surface  calculations  it  is  necessary  to  smooth  the  plan- 
form  around  the  central  and  trailing-edge-crank  sections  (Fig  1).  We  follow  the  defini¬ 
tions  in  equations  (2-7)  to  (2-9)  of  Ref  1,  apart  from  the  shift  in  origin  to  the  point  on 
the  pitching  axis  0.592c  aft  of  the  leading  apex.  It  is  necessary  also  to  fix  the  numbers 
of  spanwise  and  chordwise  terms  in  the  solutions;  again,  following  section  2,1  of  Ref  I, 
we  take  respectively  m  =  23  and  N  =  8  .  Likewise,  to  enhance  the  accuracy  of  the  solu¬ 
tions,  we  use  the  integers 

q.  =  (16,  8,  6,  4,  4,  6,  8,  16)  (2-15) 

of  Ref  7  in  evaluating  w.  from  the  integral  equation.  Nevertheless,  as  will  be  seen 
later,  the  chosen  value  of  m  is  not  large  enough  to  avoid  minor  inaccuracies. 

The  Mach  number  is  fixed  at  M  =  0.8  ,  and  it  remains  to  select  the  frequency  para¬ 
meters  v  for  which  the  complex  generalized  forces 


Q.  .(v)  =  Q!  .G)  +  ivQ'.’.(v) 

ij  ij  ij 


(2-16) 


are  required.  In  the  case  of  symmetric  sinusoidal  gusts  (j  =  5)  when  \Jj  =  0  ,  the  wave¬ 
length  in  the  streamwise  direction  is 


i 


x 


2-nc 

v 


(2-17) 


and  good  accuracy  in  superposition  for  the  step  gust  is  achieved  by  considering  the  range 
l  >  c  or,  say,  0  <  v  <  6  .  Equation  (2-17)  applies  equally  to  inclined  gusts,  for 
which  the  wavelength  in  the  spanwise  direction 


£ 


y 


2rc 

v  tan  v 


1  .572s 
v  tan  v 


(2-18) 


becomes  important.  For  example,  when  =  40°  and  v  =  6  ,  the  number  of  wavelengths 
over  the  wing  span  is  2s/£^  =  6.4  and  the  number  of  spanwise  upwash  points  m  =  23  is 
becoming  inadequate.  This  difficulty  is  intensified  at  the  larger  inclinations  i|/  ,  when 
there  is  the  further  inconvenience  that  small  values  of  v  ,  such  as  0.01,  and  fairly 
small  intervals  in  v  are  needed  to  give  adequate  definition  to  Q^,.(v)  •  For  the  Par~ 
ticular  case  =  40  ,  a  full  set  of  values  of  the  real  parts  Q!,.  is  given  in  Table  1 
for  the  36  values  of  v  used  throughout  the  present  calculations.  Except  for  the  hinge 
moment  ,  the  values  of  Q!  for  v  >  3.0  are  small  enough  in  relation  to  their 

maximum  values  to  be  almost  negligible  in  evaluating  the  Kiissner  functions,  so  that 
inaccuracies  and  wide  spacing  in  this  region  can  be  tolerated.  The  transient  character 
of  the  rolling  moment  is  typified  by  the  zero  value  of  Q^j(0)  ,  For  each  force  only  the 
first  five  stationary  values,  including  that  at  v  =  0  ,  are  thought  to  matter. 

2 . 2  Kiissner  functions  for  arbitrary  gust  inclination 

Some  detailed  attention  to  the  geometry  of  Fig  1  will  clarify  the  sequence  of  events 
after  the  wing  of  leading-edge  sweepback  A  first  meets  the  step  gust  at  time  t  =  0  . 
The  subsequent  distance  of  travel  in  mean  chords  is 


Ln 


7 


o 


(2-19) 


When  i|i  <  A  the 
the  starboard  tip 
convenient  to  use 


root  leading  edge  enters  the  gust  at  t  =  0:  when  ip  >  the  entry  of 
leading  edge  defines  t  =  0  .  For  some  purposes  in  the  analysis  it  is 


o 


0 


u(t-y 

c 


(2-20) 


where  t^  is  the  time  at  which  the  root  leading  edge  enters  the  gust.  Thus  for  \J;  >  A 
the  wing  encounters  the  gust  when 


Oq  =  0^  =  1,693  -  3.996  tan  \jj  .  (2-21) 

The  root  trailing  edge  enters  the  gust  when  a ^  =  1.690  ,  but  the  wing  does  not  become 
fully  immersed  until  the  port  tip  trailing  edge  enters  the  gust.  The  gust  penetration 
along  the  centre-line  then  defines 


Oq  =  =  2.171  +  3.996  tan  i|) 


(2-22) 


after  which  the  whole  wing  is  at  uniform  incidence  a  ,  say,  and  the  aerodynamic  loading 
2  .8 
oV  id^(x,y)  due  to  the  gust  begins  to  settle  to  a  steady  state.  The  values  of  Xq(\Jj)/c 

to  be  used  in  equation  (2-2)  for  w  =  w  are  identified  with  (o  -  0,592),  so  that 

s  e 


x0(*> 


-  0.592 

1 .  101  -  3.996  tan  ip 


(0  <  *  <  Aj) 

(A  <  ii  <  90°) 

Jc 


(2-23) 


Equations  (41a)  and  (41b)  of  Ref  2  give  alternative  expressions  for  the  time- 
dependent  force 


qiS 


•  ?lhu 


»y)  £gU,y)dxdy 


(2-24) 


integrated  over  the  planform.  In  the  present  notation  these  alternatives  are 


and 


V0) 


2a  r 

w 


Q!<-(v)  sin  \5o 


dv 


Va) 


q: 5<o)  + 


CO 

!  !  ^ 


1 


v)  cos  vo  d\i 


(2-25) 


(2-26) 


Whereas  equation  (2-25)  gives  precisely  Q^(0)  =  0  ,  equation  (2-26)  gives  an  identity 
in  the  limit  as  o  -  0  which  does  not  quite  give  zero  after  numerical  integration.  The 
present  calculations  are  based  on  equation  (2-25),  for  which  a  Fortran  program  has  been 
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written.  Because  of  the  oscillatory  behaviour  of  both  factors  in  the  numerator  of  the 
integrand  at  large  v  there  is  no  significant  error  in  truncating  the  upper  limit  of 
integration  at  a  suitably  large  value.  In  fact,  we  select  v  =  6  and  use  the  36  values 
in  Table  1  to  define  Q!^(v)  by  a  cubic  spline.  The  optimum  choice  of  upper  limit  calls 
for  judgement  as  to  where  the  inaccuracies  in  the  calculated  Q!^(v)  have  grown  to  the 
extent  that  it  is  better  to  replace  it  by  zero.  At  the  value  a  =  15  ,  say,  the  curve  of 
Q^g(a)  is  faired  into  the  asymptotic  expansion 


q;  (o)Q’  (o) 

Q.  (o)  =  Q!  r  (0)  +  — — - ¥■ -  +  0(o  J) 

lg  13  8tc“ 


(2-27) 


which  can  be  reconciled  with  equation  (38)  of  Ref  5  and  is  closely  related  to  equation 
(3-30)  of  Kef  1.  Without  this  facility  the  numerical  integration,  as  programmed,  would 
become  inaccurate  at  larger  values  of  a  because  of  the  increasingly  rapid  variations  in 


The  scope  of  the  calculations  of  Q^g(a)  includes  lift  and  the  pitching,  rolling, 
bending  and  hinge  moments  for  ten  inclinations  of  the  gust  front  in  the  range  0  <  iji  <  85° 
The  inaccuracies  associated  with  high  v  begin  to  show  themselves  for  small  o  ,  but  they 
remain  very  small  until  at  i|i  =  70°  minor  effects  obtrude  in  the  applications  of 
section  3.  Only  limited  use  can  be  made  of  the  results  for  il  =  80°  and  85°  ,  where  the 
considerations  in  equation  (2-18)  seriously  question  the  adequacy  of  m  =  23  at,  say, 
v  =  1  where  only  very  small  errors  can  be  tolerated.  The  set  of  results  for  4>  =  40° 
is  given  in  Table  2;  for  this  inclination  trivial  discrepancies  appear  as  negative  values 
of  the  bending  and  hinge  moments  on  the  trailing  (port)  side  of  the  wing  when  1  .4  <  a  < 
2.0,  and  the  hinge  moments  are  restricted  to  two  decimal  places. 

For  lift  and  for  pitching,  bending  and  hinge  moments  the  quantities  in  Table  2  are 
the  Kiissner  functions  defined  as 


K.(o)  = 


Q.  (o) 


Qig(~) 


Q^(v)  sin  vo 
Q!3  (0) v 


(i  = 1,2, 3, 5) 


(2-28) 


from  equation  (2-25)  with  truncated  integration.  By  means  of  equations  (2-14)  for  bending 
moment  and  corresponding  ones  for  hinge  moment  K^(o)  and  K^io)  are  determined  separ¬ 
ately  for  the  leading  and  trailing  sides  of  the  wing.  As  would  be  expected,  the  bending 
moment  on  the  leading  side  grows  more  rapidly  than  the  lift  or  pitching  moment.  On  the 

trailing  side  the  bending  moment  does  not  start  to  develop  until  after  o  =  -o  =  1.66 

e 

when  the  port  wing  starts  to  penetrate  the  gust.  The  rapid  growth  in  hinge  moment  on  the 
trailing  side  somewhat  anticipates  the  arrival  of  the  gust  front  at  the  aileron  corner 
at  o  =  6.05  and  is  virtually  complete  by  c  ~  =  7.18  when  the  aileron  is  fully 

immersed . 

The  rolling  moment  (i  =  4),  being  an  antisymmetric  aerodynamic  characteristic,  is 
distinctive  in  two  respects.  In  the  first  place  it  is  transient  because  Q^(0)  =  0  ; 
the  quantity  Q^g(c)  is  therefore  normalized  with  respect  to  the  steady  lift  to  give  the 
function 


8113/ 


81157 
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Q,  (o) 
Mo)  =  TT^uT 
lg 


U 

-  I  f 


Q.'.(v)  sin  vo 

4  3 _ 

q;5(o)v 


dv 


(2-29) 


The  second  distinction  is  that  the  quantity  Q'^(O)  occurring  in  equation  (2-27),  being 

the  rolling  moment  due  to  low-frequency  pitching  motion,  is  identically  zero.  As  there  is 
-2  . 

no  term  of  order  a  in  the  asymptotic  expansion,  the  steady  state  is  approached  much 
sooner  in  the  rolling-moment  column  of  Table  2. 


2 . 3  Discussion  of  results 

The  Kussner  functions,  excluding  hinge  moments,  are  used  in  idealized  active  control 
studies  at  M  =  0.8  in  section  3.  Various  aspects  of  the  aerodynamic  loading  from 
vertical  step  gusts  are  analysed  in  Figs  2  to  8,  with  particular  attention  to  the  depen¬ 
dence  on  entry  angle  p  .  Fig  2  shows  the  position  of  the  gust  front  as  the  lift  force 
reaches  a  quarter,  a  half  and  three  quarters  of  its  final  value  .  In  the  symmetric 
case  iif  =  0  the  wing  is  fully  immersed  soon  after  L  =  0.5L^  .  When  the  angle  is 
increased  to  about  20°  the  lift  reaches  0.751^  as  the  gust  front  passes  the  port  tip, 
while  for  the  most  gradual  case  ip  =  80  the  steady  state  L  =  Lq  is  approached  by  the 
time  the  whole  wing  is  immersed.  For  each  of  the  eight  inclinations  in  Fig  2  the  condi¬ 
tion  L  =  0.51.^  shortly  follows  the  instant  when  the  root  trailing  edge  enters  the  gust. 

The  function  in  equation  (2-29)  may  be  written  as 


K,  <o)  =  ~  (2-30) 

“  V 

where  £  is  the  transient  rolling  moment  due  to  the  gust.  This  quantity  is  conveniently 
plotted  in  Fig  3  against  0^  from  equations  (2-20)  and  (2-21),  and  for  each  value  of  \p 
the  maximum  occurs  fairly  close  to  c^  =  1.69  at  which  the  root  trailing  edge  enters  the 
gust.  This  maximum  grows  rapidly  as  ip  increases  from  zero.  The  steepest  gradient  in 
•£  against  occurs  soon  after  entry  when  ip  approximately  equals  the  angle  of  sweep- 

back  A.  =  23°  ,  above  which  the  curves  become  progressively  more  gentle  and  values  out¬ 
side  the  range  e  <  soon  become  negligible.  The  spanwise  centre  of  pressure  as 

a  fraction  of  wing  semi-span 

L  K4(j)? 
n  =  J.s  "  K |  (  Os 

K  (o) 

=  0.2502  t rpT  (2-31) 

K.  j  (,  O  / 


is  plotted  against  0q  in  Fig  4.  The  calculations  become  indeterminate  when  Oq  is  too 

close  to  the  entrv  value  ’  ,  and  the  curves  of  n  have  been  smoothed  to  zero  for 

e 

■p  <  .  and  to  unity  for  p  >  A,  .  The  curves  for  the  latter  range  of  all  pass  fairly 

close  to  a  point  ("Mg)  =  (0.38,0.8)  as  ~  is  decreasing  towards  zero  for  large  . 

A  few  calculations  have  been  made  to  compare  the  gust  forces  with  values  computed 
from  piston  theory^,  which  has  considerable  qualitative  merit.  The  curves  of  £/(L^c) 
from  piston  theory  in  Fig  5  and  from  lifting-surface  theory  in  Fig  3  have  the  same 
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general  trends.  Fig  6  shows  for  $  =  40°  that  piston  theory  determines  the  spanwise 
centre  of  pressure  n  well  enough.  Calculations  of  the  chordwise  centre  of  pressure 


a, 

V3 


=  0.4135 


K2(r) 

K.^) 


(2-32) 


relative  to  the  pitching  axis  by  piston  theory  correctly  predict  a  minimum  in  x/c 
against  due  to  lift  growth  near  the  leading  apex.  But  it  is  recognised  in  Figs  5 

and  6  that  piston  theory  is  a  poor  approximation  to  lifting-surface  theory  in  exaggerating 
the  peak  rolling  moments  by  factors  ranging  from  1.3  to  2.2  and  seriously  overestimating 
L/Lq  and  x/c  .  Tabulated  comparisons  of  the  lift  ratio  in  Table  3  show  maximum  discre¬ 
pancies  that  fall  from  0.43  to  0.1  I  as  *•  increases  from  0  to  70°. 

With  the  aid  of  equations  (2-14)  and  K.(c)  and  K^(a)  from  equations  (2-28)  and 
(2-29)  respectively,  the  Kiissner  functions  for  bending  moment  have  been  calculated  separ¬ 
ately  for  the  two  sides  of  the  wing.  These  are  plotted  as  B;/Bq  ant^  Bt^B0  a8ainst 
0q  in  Fig  7  for  five  values  of  ip  from  0  to  55°.  The  main  effects  are  the  steepening 
of  the  B?  curve  when  0  approaches  A.,  and  weakening  gradients  in  B  „  as  increases 
above  and  in  as  f  increases  from  zero.  The  asymptotic  expression  in  equation 

(2-27),  and  consequently  the  behaviour  for  large  ,  is  unaffected  either  by  iji  or  by 

which  side  of  the  wing  is  considered. 

Fig  8  summarizes  some  of  the  effects  of  o  and  takes  them  to  the  estimated  limit 
as  ii  tends  to  90°.  The  maximum  transient  rolling  moment  from  Fig  3,  the  spanwise  centre 
of  pressure  from  Fig  4  and  the  half-wing  bending  moments  from  Fig  7  at  the  instant  when 
L  =  0.5Lq  are  plotted  against  ip  .  The  primary  consideration  is  the  high  rate  at  which 
each  of  these  quantities  varies  while  is  small.  If  we  regard  the  range  0  <  <  20° 

as  being  statistically  less  likely  than  the  wider  range  20°  <  <  90°  ,  then  the  aero¬ 

dynamic  characteristics  for  this  upper  range  are  the  more  representative.  They  will  be 
modified  by  gust  length,  as  we  shall  find  in  section  4,  but  the  maximum  rolling  moment 
will  remain  substantial. 

3  NEUTRALIZATION  OF  GUST  FORCES 

In  general  the  response  to  asymmetric  gust  entry  introduces  all  the  lateral  degrees 
of  freedom.  In  the  present  studies  the  ailerons  are  used  to  neutralize  the  rolling 
moment  due  to  the  gust.  This  idealization  leads  to  the  suppression  of  motion  in  a  single 
antisymmetric  degree  of  freedom  in  roll,  and  the  question  is  how  the  asymmetry  is  likely 
to  influence  control  demand  or  structural  loading  in  the  context  of  active  control.  As 
in  Ref  1,  perfect  knowledge  of  the  aerodynamic  input  from  the  gust  is  assumed  without  the 
practical  consideration  of  sensors  or  control  laws.  We  can  devise  two  numerical  studies 
in  search  of  any  adverse  effects  of  the  gust  inclination  . 

The  first  example  is  an  idealized  study  of  load  alleviation  with  complete  disregard 
of  the  modes  of  heave  and  pitch.  The  simplified  problem  is  to  determine  separate  move¬ 
ments  of  the  two  ailerons  so  as  to  eliminate  si mul taneously  the  aerodynamic  root  bending 
and  rolling  moments  throughout  the  gust  encounter.  The  influence  of  ^  on  control 
demand  can  then  be  analysed. 


Cr 
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The  second  numerical  study  simulates  rise  control  for  a  rigid  wing  with  the  suppres¬ 
sion  of  the  forces  which  would  cause  motion  in  heave,  pitch  and  roll  throughout  the  gust 
encounter.  To  cater  for  the  three  degrees  of  freedom,  one  additional  control  surface  is 
required.  This  third  control  is  conceived  as  a  tailplane  or  canard,  which  is  specified 
simply  as  a  lift  force  at  a  prescribed  location  on  the  centre-line  x  =  or  x  =  x  . 
The  aerodynamics  of  this  surface  are  ignored  in  three  respects,  wing  interference,  its 
gust  loading  and  the  lag  in  growth  of  lift  as  it  is  deployed.  With  this  simplification 
the  iterative  scheme  of  calculation  embodied  in  equations  (5-5)  to  (5-10)  of  Ref  1  will 
suffice  to  determine  the  time  history  of  the  two  aileron  angles  and  the  control  lift 
or  .  The  objective  in  section  3.3  is  to  analyse  these  quantities  and  the  residual 

bending  moment,  which  plays  the  role  of  structural  penalty. 

3 . I  Wagner  functions  for  control  surfaces 

The  dependence  of  time-dependent  aerodynamic  forces  on  the  complete  history  of  wing 
and  control-surface  motion  is  treated  according  to  linearized  theory  in  Ref  1 .  Several 
kinds  of  hereditary  function  are  discussed  in  that  paper,  and  in  the  present  application 
we  shall  identify  the  'hereditary  velocity  factor’  of  Ref  1  as  a  generalization  of  the 
Wagner  function  of  Ref  3.  Equation  (2-22)  of  Ref  1  is  therefore  rewritten  as 


Wij(c) 


2 


QV.(v) 


QV . (a>) 


Q!.(0) 

1J 


cos  vo  dv 


There  is  an  alternative  expression  similar  to  that  for  the  Kiissner  function  in  equation 
(2-28),  but  equation  (3-1)  is  preferable  for  the  purposes  of  calculation.  Section  3.3 
and  Appendix  B  of  Ref  1  describe  how  this  is  evaluated  by  subdividing  the  infinite  range 
of  integration  into  three  parts. 

With  application  to  active  controls  the  upward  displacement  in  the  mode  j  is 
given  an  arbitrary  smooth  time  dependence 


z  =  q^(T)Zj(x,y) 


(j  =  1,2,3) 


where  t  =  Ut/c  and  z.  corresponds  to  equations  (2-8).  We  shall  be  concerned  with  the 
symmetric  aileron  motion  j  =  3  and  antisymmetric  aileron  motion  j  =  4  corresponding 
to  w^  in  the  second  of  equations  (2-10),  To  avoid  inconsistency  in  notation,  the 
equivalent  of  equation  (3-2)  for  j  =  i  is  written  as 


VT)(lf])z3(x’y)  • 


Then,  following  equation  (2-21)  of  Ref  I,  we  shall  use  W.j(a)  of  equation  (3-1)  to 
obtain  contributions  to  the  aerodynamic  force  in  mode  i 


Q  •  (  t)  =  QV • (”) 
i  i  J 


dq. (t) 


T 


dq.(tn) 


wij(T"To)dTo 


i 


12 


The  Wagner  functions  for  lift  and  the  pitching,  hinge,  rolling  and  bending  moments 
(i  =  1  to  5  respectively)  have  been  calculated  for  aileron  motion  of  the  configuration  in 
Fig  1  at  M  =  0.8  and  are  reproduced  in  Table  4.  As  regards  lift,  the  trailing-edge  con¬ 
trol  has  a  more  sluggish  aerodynamic  response  than  the  other  types,  leading-edge  control 
and  all-moving  tip,  considered  in  Ref  1.  The  initial  value  Wj3(0)  =  0.343  is  well  below 
the  value  0.5  associated  with  the  original  Wagner  function  for  heaving  velocity  in  two- 
dimensional  incompressible  flow.  But  the  asymptotic  expansion  for  large  a  in  three 
dimensions 


W..(a) 


Ql.i(0)Qi2(0) 

8ttQ^(0)o2 


0(a“3) 


(3-5) 


gives  a  more  rapid  approach  to  the  steady  state  than  in  two  dimensions  when  there  is  a 
term  of  order  a  '  .  Thus  the  value  W^(o)  =  0.978  when  a  =  10  compares  with  the 
value  0.937  of  the  original  Wagner  function.  The  pitching  moment  for  the  particular  axis 
responds  much  more  rapidly  than  lift  and  gives  Woj(10)  =  0.992  .  The  bending  moment  has 
intermediate  characteristics  so  that  for  all  a 


W13(a)  <  W53<a)  <  W,3(a)  .  (3-6) 

The  respective  Wagner  functions  W^(o)  and  W33(a)  for  the  rolling  moment  due  to 

antisymmetric  ailerons  and  for  the  symmetric  hinge  moment,  approach  the  steady  state  much 

more  rapidly  for  different  reasons.  Like  the  Kussner  function  K^(o),  W^(o)  has  no  term 

of  order  o  2  in  the  asymptotic  expansion  (3-5)  because  Q'  (0)  =  0  ;  it  therefore 
-3  ^ 

behaves  like  1  +  0(o  )  in  Table  4.  The  clue  to  the  behaviour  of  the  symmetric  hinge 

moment  is  found  in  the  values  of  -Q'33(v)  in  Table  2  of  Ref  1,  which  are  almost  constant 
over  the  range  0  <  v  <  2  .  The  factor  £q'33(v)  -  Q'^i00)]  in  the  numerator  of  equation 
(3-1)  is  not  especially  small,  roughly  25%  of  Q33(0)  ,  but  the  constancy  of  that  factor 
in  association  with  cos  vo  reduces  the  magnitude  of  the  contribution  of  the  integral 
to  0.03  or  less  for  a  >  1  .  With  this  feature  and  its  large  initial  value  W33(0)  = 
1.908,  the  hinge-moment  function  is  strikingly  different  from  the  original  Wagner 
function  and  is  included  in  Table  4  for  its  intrinsic  interest.  The  other  four  are 
required  subsequently. 

3 .2  Load  alleviation  and  control  demand 

This  exercise  requires  the  Kussner  functions  K, (o)  and  K^Co)  for  the  whole 
range  of  and  the  Wagner  functions  W^(o)  and  W33(c)  .  The  objective  is  to  deter¬ 
mine  the  aileron  motions  q„(x)  and  q,(x)  in  equations  (3-2)  and  (3-3)  so  as  to  give 

J  *4 

quantities  Q^(x)  in  equation  (3-4)  to  neutralize  the  rolling  and  bending  moments  due  to 
the  step  gust.  From  the  definitions  in  equations  (2-28)  and  (2-29)  the  integral  equations 


V")K5(t) 


Q'5'3(») 


dq3(x) 


dr 


+  Q^j(0) 


dci3(V 

dx  ^ 


W53(T_T0)dT0 


=  0 


(3-7) 
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dq, (t) 

Qlg(“)VT)  +  Q44(M)  -Tr~  + 


r  dW 

Q44<°>  J  W, 


44(T_T0)dT0  =  ° 


have  to  be  solved  for  q^Ct)  amd  q4(x)  .  The  two  equations,  being  symmetric  and  anti¬ 
symmetric  respectively,  are  independent  and  can  each  be  solved  by  the  iterative  process 
based  on  equation  (5-5)  of  Ref  1.  Equation  (3-7)  is  written  in  the  form  of  a  linear 
differential  equation  for  q^n  *^(x) 


Q»<->  v->y> 

^  (  }  Q'  (0)  dr  Q'  (0)  J  dt  \W53(t  t0)  )dT0  ’ 


Q’53(0)  dt 


q;3(°) 


where  q^  (r)  is  a  previous  approximation  to  q  (x)  .  With  initial  values  ql^(x)  =  0 

J  d  (9)  d 

the  process  is  found  to  converge  satisfactorily  to  q3(i)  *  q^  (x)  after  nine  itera¬ 
tions.  Equation  (3-8)  is  solved  for  q^(x)  similarly. 

The  following  values  are  used  in  the  calculations: 


q  («)  =  -  a  Q',(0)  =  43.17a  for  all  4) 

5g  g  55  g 


Q’3(0)  =  -  5.386  ,  Qj3 (“)  =  -  0.178 


(3-10) 


Q)g(») 


a  Q'  (0)  =  25.37a  for  all  4 

g  15  g 


q;4(°)  =  -5.113,  Q44(”)  =  -0.178  J 

where  ag  is  the  upwash  angle  in  radians  associated  with  the  gust.  For  the  steady  state 
at  large  time  equation  (3-7)  or  (3-9)  gives 


Q5k(“> 

q3(<”)  ~  ~  Q’jiO)  =  8-015ag  * 


(3-11) 


and  the  corresponding  limit  from  equation  (3-8)  is  q,  (»)  =  0  because  K^i")  =  0  as  in 
Table  2.  We  write 


d^(T)  =  ~  q3(t)  -  q4(x) 


^ t ( -t )  =  -  q3(x)  +  q4(x) 


(3-12) 


for  the  downward  (nose-up)  deflections  of  the  ailerons  in  radians  on  the  leading  and 
trailing  sides  of  the  wing.  The  limiting  conditions 


S ,  (“)  =  6  ( »)  =  -  8.0)5a 

1  t  g 


(3-13) 


14 


imply  nose-down  deflections  of  both  ailerons  to  counteract  the  upward  bending  moment  due 
to  an  up-gust  represented  by  positive  a . 

From  the  calculations  of  6{(r)  and  <5 1 <[ r )  for  various  gust  inclinations  in  the 

range  0  <  \Ji  <  80°  ,  a  few  illustrative  curves  of  -6  /a  and  -<5  /a  against  the  par a- 

t  g  t  g 

meter  of  equation  (2-20) 


t  _^0 
c 


(3-14) 


are  drawn  in  Fig  9.  The  asymptotic  from  equation  (3-13)  is  approached  by  the  time  =10 

or,  if  >  55°  ,  as  soon  as  the  wing  is  fully  immersed  in  the  gust.  The  curves  have 

positive  slopes  everywhere,  and  the  main  influence  of  4 i  is  to  steepen  the  curve  on  the 
leading  side  when  i|i  =  25°  ,  which  is  close  to  the  condition  tp  =  A  when  the  whole  lead¬ 
ing  edge  of  the  leading  side  enters  the  gust  at  once.  Above  this  on  the  leading  side  and 
for  all  tp  on  the  trailing  side  the  gradients  slacken  with  increasing  4/  ,  and  there  are 
no  adverse  effects  of  asymmetric  entry  into  a  step  gust.  The  primary  reasons  are  the 
increasing  distance  of  travel  from  gust  entry  to  full  immersion  as  increases,  (a^-o  ) 
mean  chords  from  equations  (2-21)  and  (2-22),  and  the  rapid  response  in  rolling  moment 

from  antisymmetric  aileron  movement,  which  is  characterized  by  the  disappearance  of  the 

-2 

a  term  in  the  asymptotic  expansion  (3-5)  for  W^(a)  . 

3 . 3  Ride  control  and  residual  bending  moment 

In  Ref  I  there  are  calculations  of  symmetric  aileron  motion  to  neutralize  the  lift 
due  to  symmetric  entry  into  a  step  gust.  For  the  present  wing  of  Fig  1  at  M  =  0.8  ,  it 
is  concluded  in  Ref  1  that  the  aileron  motion  6(t)  stabilizes  about  three  times  more 
quickly  than  the  growth  of  lift  Kj (c)  .  The  objective  now  is  to  extend  that  calculation 
to  three  degrees  of  freedom  in  heave,  pitch  and  roll  and  to  study  the  consequences  for 
wing  bending  moment.  We  use  the  Kussner  functions  (o)  ,  ^(c),  K^(o)  and  K^(o)  from 
Table  2  for  iJj  =  40°  and  the  corresponding  functions  for  other  gust  inclinations  and  the 

Wagner  functions  Wjj(o),  W^^(a) ,  W^(c)  and  W  ^(o)  from  Table  4. 

The  control  motions  are  split  into  symmetric  and  antisymmetric  aileron  deflections 
q^(r)  and  q^(t)  ,  as  in  section  3.2;  for  the  first  example  we  introduce  a  tailplane  to 
provide  auxiliary  lift  L^,(t)  at  a  supposedly  fixed  location  x  =  x^  .  These  three 
unknowns  and  the  residual  bending  moment  B(t)  can  be  obtained  from  separate  calcula¬ 
tions.  The  antisymmetric  term  q^(x)  satisfies  equation  (3-8)  and  is  solved  iteratively 

from  the  linear  differential  equation 


(n+ 1 ) ,  ,  .  <$4<-> 
q4  (T)  Q^(0)  dt 


q;4(o> 


t  ,  (n)  ,  . 


44(t_T0) 


'ho 


(3-15) 


precisely  as  in  section  3.2.  To  obtain  q^it)  ,  the  corresponding  process  is  used  for  the 
pitching  moment  about  the  axis  x  =  x^.  ,  to  which  the  tailplane  does  not  contribute.  The 
appropriate  Kussner  and  Wagner  functions  K^.(o)  and  WTj(o)  are  simply  derived  from  the 
definitions  in  equations  (2-28)  and  (3-1).  Remembering  that  for  i  =  I  and  2 
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Qig(“)  =  -^“g  - 

we  obtain 

QJ5(0)K2(o)  -  (xt/H)Q;5(0)K,(o) 

Vo)  =  Q^5(0)  -  (xt/C)Q;5<0) 

and 

Q23(0)W23(o)  "  (xt/?)Q!3(0)WI3(o) 
wt3(o)  =  q^3(o>  -  (xt/oq;3(o) 


(3-16) 


(3-17) 


(3-18) 


The  iterative  solution  for  q3(t)  follows  from  the  equation 


q3  (t)  Q^3(0)  di 


qt8(co)kt(t)  [ 

Q^(°)  -J 


T  dq<n)(TA> 


dt. 


|wt3(t-t0)  I  }dT0 


(3-19) 


where  from  Table  1  of  this  paper  and  Table  2  of  Ref  1 


Q^3<«) 

"  Q23(~>  - 

c 

Q'.'3(-) 

XT 

-  0.051  — 

c 

-  0.073 

QTg(lo) 

-  «25(0)  " 

q;5(0) 

XT 

=  25.37  -4- 

c 

-  10.49 

a 

S 

c 

Q’3(0) 

■  <^3(0)  " 

^T 

c 

q;3(0) 

XT 

=  1.788— 

c 

-  1.884 

Having  obtained  q3<t)  and  q^(t)  after  nine  iterations,  we  evaluate  Lt(t)  and 
B(x)  from  equations  (2-28)  and  (3-4).  For  the  lift,  which  must  be  equated  to  zero. 


Lt(t) 


a  O'  (O)K.(t)  +  Q"  («) 


dq3(t) 


puV  ^'5 . ' 


13  dt 


f  dq3(T0) 

+  -i~W.3^-VdT0  "  °  (3-2l) 

6  0 


and  for  the  residual  bending  moment 


B(T) 

T,2-3 

pU  c 


-  agQ’5(0)K5(T)  +  q^(-)  dT 


dtJ3(T)  .  At  _  j  dq3(V 


+  Q33 (0) 


dt_ 


W53(x-T0)dT0  .  (3-22) 


Because  the  rolling  moment  is  neutralized,  the  bending  moment  consists  of  equal  contribu¬ 
tions  B  and  B^  from  the  two  sides  of  the  wing.  The  individual  aileron  angles 
and  5  are  given  by  equations  (3-12). 

A  epical  set  of  results  for  \|i  =  40°  and  x^,  =  3c  is  listed  in  Table  5  for  a 

unit  down-gust  a  =  -1  .  The  control  angles  6.  and  6  are  therefore  positive  to 

g  *  _  c 

neutralize  the  negative  lift  due  to  the  gust.  The  tailplane  lift  is  of  the  opposite  sign, 
because  the  centre  of  lift  due  to  the  gust,  x/c  in  Fig  6,  lies  forward  of  that  due  to 
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the  aileron  motion.  The  bending  moment  is  very  roughly  proportional  to  (6  + £  ) .  The 
normalized  values  of  lift  and  bending  moment  are  L^/L^  anc*  B/Bq  w^ere  c^e  positive 
quantities 


lq  =  ^22\q;5(°) 


(3-23) 


Bfl  =  pU2S3agQ'5(0) 


are  the  steady  forces  on  the  wing  at  unit  incidence.  The  limiting  values  as  a  tends  to 
infinity  follow  from  equations  (3-12),  (3-15)  and  (3-19)  to  (3-23)  which  give 


QT g(”) 

"  q3(“)  "  Qt3(0) 


[25.37(5^/5)  -  10.49V 


1.788(xt/c)  -  1.884 


(3-24) 


An  interesting  feature  of  Table  5  is  that  each  of  the  quantities  6^,  S^,  -L^/Lq 
and  B/B0  reaches  a  maximum  near  a  -  7  before  falling  to  its  steady  limit.  This  is 
associated  with  the  lag  in  the  growth  of  lift  when  a  trailing-edge  control  is  deployed 
rapidly.  The  overshoot  in  aileron  angle  arises  because  the  high  limiting  value  in 
equation  (3-24)  when  x^,  =  3c  demands  a  control  rate  high  enough  to  need  a  significant 
extra  boost  to  cope  with  the  consequent  lag  in  lift. 

To  change  the  tailplane  to  a  canard  it  is  sufficient  to  replace  the  positive 
quantity  stj,  and  the  subscript  T  ,  wherever  they  have  appeared,  by  a  negative  quantity 
xc  and  the  subscript  C  .  The  introduction  of  the  extra  symbol  is  tc  emphasize  that 
there  is  an  intermediate  regime  where  singularities  would  inhibit  the  consideration  of 
active  control  and,  in  particular,  values  where  equations  (3-17),  (3-18)  and  (3-24)  to 
(3-26)  would  fail.  Corresponding  calculations  for  ij;  =  40°  and  xc  =  -3c  in  Table  6 
have  features  in  contrast  to  those  in  Table  5  with  the  tailplane  as  third  control. 

Because  the  lift  forces  from  the  canard  and  the  ailerons  now  act  on  opposite  sides  of  that 
from  the  gust,  at  x/c  in  Fig  6,  the  pitching  moment  can  be  neutralized  with  smaller 
individual  lift  contributions  of  the  same  sign.  In  consequence  the  bending  moment  is 
substantially  reduced  and  the  active  canard  is  seen  to  be  favourable  to  the  structure. 

This  is  apparent  from  the  static  considerations  inherent  in  equations  (3-24)  to  (3-26) . 

The  smaller  demand  for  control  rate  removes  the  overshoot  found  in  Table  5  with  the  tail- 
plane,  so  that  with  the  canard  each  of  the  quantities  5p,  6^,  and  B/Bq 

within  1%  of  its  final  value  as  soon  as  o  >  9  . 
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The  calculations  covering  variations  in  i,  x^,  and  x^  are  summarized  in  Figs  10 
ti  !5.  Fig  10  gives  the  behaviour  of  ‘  =  f  =  -  ^  and  or  Lc  in  symmetric  cases 

that  differ  from  those  in  Tables  5  and  b  in  that  0=0  instead  of  tp  =  40°  .  The  limit¬ 
ing  values  in  equations  ( 3—24)  and  (3-25)  are  unchanged,  but  the  overshoots  in  the  tail- 
plane  case  x^  =  3c  are  somewhat  more  marked.  The  canard  case  x^,  =  -3c  is  special, 
because  the  steady  state  is  virtually  reached  when  =  6  .  The  sluggishness  in  growth 
of  lift  due  to  aileron  motion  requires  additional  control  rate  that  is  just  insufficient 
to  give  any  overshoot.  Fig  11  gives  further  information  on  the  effect  of  iJj  on  control 
demand  in  the  tailplane  case  Xj,  =  3c  .  The  most  rapid  control  rate  appears  at  p  =  25°  , 
being  close  to  y  =  V.  •  but  the  increase  is  not  large  enough  to  aggravate  the  overshoot. 

A  small  overshoot  remains  on  both  sides  of  the  wing  at  ip  =  55°  ,  At  ip  -  70°  the 
inaccuracies  in  the  basic  lifting-surface  calculations,  mentioned  below  equation  (2-18), 
are  thought  to  be  responsible  for  the  undulations  in  and  <5^  ,  which  are  only 

apparent  for  the  first  few  mean  chords  of  travel  after  gust  entry  at  o^  =  -9.3  . 


Fig  12  and  the  upper  diagram  of  Fig  13  concern  the  influence  of  the  location  of  the 

third  control,  which  has  been  calculated  for  p  =  0  and  40°  .  Curves  of  -6  /a  and 

1  S 

-<5t/a  against  an  for  the  asymmetric  case  are  given  in  Fig  12,  and  they  virtually 


'0 

coalesce  near  Cq  =  6  in  the  three  cases  x^,  =  -2c  and  x^  =  2c  and  4c  .  With  the 

shorter  tail  arm  the  overshoot  in  -5  ,/a  is  9%,  which  falls  to  2%  when  x„  =  4c  .  In 

*  S  ^ 

the  absence  of  overshoot  for  the  canard  control  at  x^,  =  -2c  the  maximum  control  angles 

equal  the  limiting  one  -^(“O/a  =  11.2  from  equation  (3-24),  which  provides  the  curve 

of  -6  /a  in  Fig  13  with  canard  lift  whatever  the  value  of  p  .  With  tailplane  lift 
max  g 

there  is  slight  dependence  of  "®max/a  on  '!’  because  the  overshoot  is  less  at  =  40 
than  at  =  0  .  More  significant  by  far  is  the  reduction  by  a  factor  of  2,  say,  when 
the  tailplane  is  replaced  by  a  canard. 


If  S  is  taken  as  a  measure  of  control  demand,  the  influence  of  asymmetric  gust 

max 

entry  is  marginally  beneficial  in  the  tailplane  case.  The  extent  of  this  small  allevia¬ 
tion  for  x_  =  3c  is  shown  in  the  lower  diagram  of  Fig  13,  where  -6  /a  falls  from 
T  6  b  max  g 

19.9  to  the  final  value  18.8  as  if  increases  from  0  to  90  .  If  the  residual  root  wing 

bending  moment  is  regarded  as  a  structural  penalty  in  achieving  ride  control,  then  Fig  14 

again  shows  the  influence  of  asymmetric  gust  entry  to  be  marginally  beneficial;  whether  a 

tailplane  or  a  canard  is  used,  B  falls  a  little  as  increases  from  zero.  The 

max 

overriding  consideration  for  all  p  is  that,  with  canard  lift  as  opposed  to  tailplane 
lift,  the  residual  bending  moment  is  reduced  by  a  factor  of  about  3J.  If,  however,  the 
maximum  control  rate 

{max  =  [^("^)]  (3_27 

max 


is  regarded  as  the  primary  influence  of  the  gust  inclination,  then  the  peak  values  near 
„  =  A(  =  23.0°  for  the  aileron  on  the  leading  side  of  the  wing  appear  as  the  single 
adverse  effect  of  asymmetric  gust  entry  to  emerge  from  the  active  control  studies  for 
step  gusts.  Fig  15  shows  the  extent  of  this  phenomenon  both  for  ride  control,  achieved 
with  active  ailerons  and  tailplane  lift,  and  for  load  alleviation  as  discussed  in 
section  3.2. 
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4  INFLUENCE  OF  GUST  LENGTH 

Effects  of  gust  length  are  considered  in  relation  to  (1-cosine)  profiles,  which 
correspond  to  a  half  cycle  of  a  sinusoidal  gust.  If  x^  denotes  the  streamwise  length 
of  the  gust,  then  in  place  of  equation  (2-2)  we  have  an  instantaneous  upwash  angle 


w(x,y) 

U 


=  -£<1  -  cos  — (ut  -  x  +  y  tani^  +  x-(i)j) 


Ut  -  x  +  y  tan  if'  +  xq(0)  <  0 


0<Ut-x+y  tan  i*/  +  Xq(^)  <  x^ 
Ut-x+y  taniJj  +  Xg(^)  >  x^ 


As  x^  tends  to  zero,  the  middle  region  disappears  to  give  a  step  gust  of  incidence 
met  by  the  wing  at  non-dimensional  time  t  =  0  .  The  (1-cosine)  gust  consists  of 
elementary  gusts  of  incidence 


na 

6a  =  si 

2o 


x 

c  =  -& 

8  c 


‘"Q-o 


at  times  0  <  T_  <  o 
0  g 


By  superposition,  therefore,  any  function  of  o  corresponding  to  a  step  gust  of  incidence 
a  and  a  particular  inclination  'p  ,  say  the  rolling  moment  X(o)  ,  can  be  generalized  to 


•C(c',og)  = 


j  lir  sin(-r)£{a  ~  To)dTo 

0  8  \  8  / 


where  a  is  a  or  o  ,  whichever  is  the  smaller, 
m  g 

By  means  of  equation  (4-4)  it  is  a  simple  calculation  to  extend  the  results  of  Fig  3 
for  X/LgC  in  an  inclined  step  gust  to  corresponding  functions  for  (1-cosine)  gusts  of 
different  length.  For  the  particular  gust  inclination  &  =  40°  at  M  =  0.8  ,  Fig  16 
shows  eight  curves  of  I/L^c  against  o  for  gust  lengths  in  the  range  0  <  o  <  60  . 

The  inset  diagram  of 


t  ■  *['  '  “©] 


0  <  1 

a 

g 


in  Fig  16  indicates  the  gust  incidence  at  the  starboard  tip  leading  edge  from  the  instant 

o=0  when  it  first  enters  the  gust.  The  maximum  rolling  moment  occurs  after  travel  of 

about  3.5c  beyond  the  point  where  a  =  and  o  =  ;  the  maximum  value  decreases 

slowly  at  first  as  o  increases,  but  when  a  =20  it  has  fallen  to  34%  of  its  value 
8  g  o 

for  the  step  gust.  Fig  17  includes  the  curve  of  X  /L„c  against  a  for  =  40  , 

_  j  IHclX  0  g 

which  is  roughly  proportional  to  for  large  gust  lengths  and  remains  significant 
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up  to  c  =  60  .  Similar  calculations  of  lift  L(o,o  )  for  ib  =  AO  and  0  <  o  <35 
8  8  8 
are  listed  in  Table  7,  where  the  greatest  spread  of  values  occurs  at  a  =  10  .  The  values 

of  a  at  which  L/L„  =  0.25,  0.50  and  0.75  are  obtained  for  0  <  o  <60  and  are  plotted 
0  g 

as  nearly  straight  lines  in  Fig  17.  The  corresponding  spanwise  centres  of  pressure  n 
for  constant  L/L^  and  the  maximum  rolling  moment  have  the  same  general  behaviour  against 
0^  in  Fig  17.  Their  behaviour  against  i |/  is  similar  to  that  shown  in  Fig  8, 

The  graphs  in  Figs  18  to  °0  illustrate  the  softening  influence  of  o^  on  the  results 
of  the  active  control  studies.  The  control  angles  and  tailplane  lift  in  Figs  18  and  19  are 
obtained  as  functions  of  o  and  o^  from  integrals  analogous  to  equation  (4-4).  Fig  18 
shows  the  gentler  gradients  in  the  two  aileron  angles  when  the  load-alleviation  exercise 


of  section  3.2  is  repeated  for  tji  =  40  and  increasing  values  of 


from  0  to  35. 


Fig  19  illustrates  the  effect  of  on  the  ride-control  exercise  of  section  3.3  in  the 

symmetric  case  =  0  with  the  tailplane  at  x^,  =  3c  .  As  would  be  expected,  the  over¬ 

shoot  of  about  6%  in  -S/a  and  -Lt/Lq  in  the  case  of  a  step  gust  decreases  as 
increases  and  has  almost  disappeared  when  o  =  20  .  The  complete  asymmetric  active  con- 

8  Q 

trol  solution,  including  the  bending  moment,  is  given  for  M  =  0.8,  =  10  and  ip  =  40 

in  Table  8.  A  slight  overshoot  in  all  quantities  is  found  near  o  =  16  .  The  influence 

of  a  can  be  studied  by  comparing  Table  5  for  o  =0  with  Table  8  for  o  =  10  . 

8  8  8 
When  a  >  25  the  differences  no  longer  exceed  one  in  the  last  decimal  place.  Fig  20 


gives  various  curves  of  maximum  control  rate  from  equation  (3-27)  against 


Although 


-1 


these  show  the  usual  behaviour  proportional  to  a  '  for  large  a  ,  the  percentage 
.  8  8 
decreases  in  <5^^  in  the  range  0  <  a  <20  are  greater  than  the  corresponding  decrease 

in  maximum  transient  rolling  moment  in  Fig  17. 


It  is  appropriate  to  .  onsider  disturbances  like  small  step  gusts  in  discussing  the 
implications  of  asymmetric  gust  entry  for  problems  of  active  control.  But  the  maximum 
transient  rolling  moment  needs  to  be  evaluated  with  reference  to  a  typical  (1-cosine) 
gust  in  the  airworthiness  requirements.  This  aspect  of  asymmetric  gust  entry  is  con¬ 
sidered  as  part  of  a  broader  discussion  in  section  5. 


5  SIGNIFICANCE  OF  ASYMMETRIC  GUST  ENTRY 

Asymmetry  in  gusts  poses  a  variety  of  problems.  The  present  studies  only  concern 
horizontal  surfaces  in  vertical  gusts,  but  the  airworthiness  requirements  recognise  the 
corresponding  case  of  a  fin  in  a  horizontal  gust.  Specifications  in  terms  of  a  (1-cosine) 
gust  of  velocity  66  ft/s  (20  m/s)  and  length  25c  have  evolved  empirically.  This  particu¬ 
lar  (l-cosine)  gust  is  taken  to  be  representative  of  the  spectrum  of  current  airworthiness 
requirements.  More  severe  velocities  can  occur  in  thunderstorms  and  jet  streams,  and  it 
is  thought  that  shorter  gradient  distances  occur  naturally  in  vortex  flows  induced  in  the 
atmosphere.  Asymnetric  spanwise  loading  may  arise  from  lateral  variations  in  gust 
velocity,  especially  in  relation  to  natural  vortex  flows  and  aircraft  wakes.  The  parti¬ 
cular  aspect  of  asymmetric  gust  entry  is  one  such  problem  where  some  clarification  is  now 
possible. 

The  piesent  numerical  studies  have  developed  in  the  context  of  active  controls  and 
fairly  small  vertical  gusts.  The  two  idealized  applications  are  to  load  alleviation, 
where  the  ailerons  are  used  independently  to  eliminate  the  rolling  and  root-bending 
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moments,  and  to  ride  control,  where  lift  and  the  pitching  and  rolling  moments  are  neutra¬ 
lized  by  combined  motion  of  the  ailerons  and  a  lift-producing  tailplane  or  canard.  The 
influence  of  the  direction  of  entry  into  the  gust  is  analysed,  and  the  maximum  control 
angles  and  residual  loads  are  little  changed.  In  the  case  of  ride  control  with  a  tail- 
plane,  an  overshoot  in  control  angle  is  demanded  (Fig  10)  and  there  is  a  corresponding 
peak  in  the  structural  penalty  measured  by  the  bending  moment  (Fig  14).  Both  these  over¬ 
shoot  effects  diminish  as  the  direction  of  entry  becomes  more  inclined  (Figs  13  and  14) 
and  as  the  gust  length  increases  (Fig  19);  asymmetric  gust  entry  is  marginally  beneficial 
in  this  respect.  The  only  adverse  effect  is  in  relation  to  maximum  control  rate  demand 
(Fig  15),  which  rises  to  a  peak  for  the  swept  wing  when  the  direction  of  entry  is  normal 
to  the  leading  edge.  More  precise  numerical  studies  could  be  made  by  considering  the 
time-dependent  aerodynamics  of  lifting  surfaces  in  tandem.  But  the  calculated  influence 
of  asymmetric  gust  entry  is  not  significant  enough  to  warrant  an  extended  investigation 
into  interfering  surfaces  and  the  Kiissner  effect  on  the  tailplane  or  canard. 

If  the  present  results  are  applied  to  the  (1 -cosine)  gust  of  velocity  66  ft/s 
(20  m/s),  it  is  easily  verified  that  the  calculated  demands  on  control  angle  and  control 
rate  are  impractically  high.  Even  for  the  (1-cosine)  gust  of  length  100  ft  (30.5  m)  ,  the 
maximum  control  rate  in  Fig  20  with  neutralized  wing  forces  and  tailplane  lift  would 
exceed  1000°/s  and,  whatever  the  gust  length,  the  aileron  angles  in  Fig  19  would  be 
approaching  90°.  It  is  more  to  the  point  that  asymmetric  gust  entry  gives  rise  to  a  sub¬ 
stantial  rolling  moment,  whose  peak  value  grows  rapidly  as  the  entry  angle  increases  from 
zero  (Figs  3  and  8).  For  a  step  gust  of  inclination  =  40°  the  maximum  rolling  moment 
is  already  86%  of  its  limit  at  tp  =  90°  with  precisely  half  the  wing  immersed  in  the 
gust.  A  further  reduction  from  86%  to  25%  is  found  when  the  gust  length  is  increased  to 
25c  (Figs  16  and  17).  Nevertheless,  the  maximum  transient  rolling  moment  due  to  a 
(l-cosine)  gust  of  this  length  and  velocity  66  ft/s  (20  m/s)  is  equivalent  to  that  pro¬ 
duced  by  static  aileron  deflections  of  +3.4°  at  M  =  0.8  ,  or  larger  deflections  roughly 
proportional  to  I/M  at  lower  flight  speeds. 

The  present  evidence  suggests  that  the  traditional  aircraft  response  calculations 
with  symmetric  gust  entry  cater  for  the  severest  structural  loads.  In  practice  there  will 
be  spanwise  asymmetries  as  well,  either  from  the  nature  of  the  gust  itself  or  on  account 
of  the  direction  in  which  it  is  approached.  The  control  rate  demands  can  then  be  greater, 
and  the  acquisition  of  bending  moments  will  differ  on  the  port  and  starboard  wings.  The 
problem  of  aircraft  control  in  gusts  is  oversimplified  if  all  gusts  are  assumed  to  be 
symmetrical . 

6  CONCLUSIONS 

(1)  Lifting-surface  calculations  show  that  the  transient  rolling  moment  due  to  asymmet¬ 
ric  gust  entry  builds  up  rapidly  as  a  function  of  the  inclination  of  flight  direction  to 
the  normal  to  the  gust  front. 

(2)  From  particular  calculations  at  M  =  0.8  the  maximum  rolling  moment  experienced  in 
a  vertical  (l-cosine)  gust  of  velocity  66  ft/s  (20  m/s)  and  length  25c  is  large  enough  to 
be  equivalent  to  that  produced  by  static  aileron  deflections  of  i3.4°. 
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Table  I 


OSCILLATORY  DATA  FOR  SINUSOIDAL  GUSTS  OF  VARYING  WAVELENGTH 
AND  OF  INCLINATION  'ji  =  40°  at  M  =  0 . 8 


V 

Lift 

Pitching 

moment 

Rolling 

moment 

Bending 

moment 

Hinge 

moment 

-Q’l5 

~Q25 

-Q55 

~Q35 

0.000 

25.37 

10.49 

0.00 

43.17 

0.0082 

0.010 

25.32 

10.47 

0.02 

43.09 

0.0082 

0.020 

25.19 

10.42 

0.09 

42.84 

0.0082 

0.035 

24.84 

10.27 

0.27 

42.21 

0.0081 

0.050 

24.33 

10.06 

0.54 

41.28 

0.0079 

0.075 

23.17 

9.57 

1.21 

39.18 

0.0076 

0.100 

21  .71 

8.96 

2.12 

36.51 

0.0072 

0.125 

20.03 

8.24 

3.24 

33.43 

0.0067 

0.150 

18.19 

7.46 

4.55 

30.07 

0.0061 

0.180 

15.88 

6.47 

6.31 

25.83 

0.0054 

0.210 

13.55 

5.45 

8.22 

21.55 

0.0047 

0.250 

10.51 

4.12 

10.87 

15.98 

0.0037 

0.300 

7.00 

2.56 

14.10 

9.60 

0.0025 

0.350 

3.93 

1  .20 

16.93 

4.13 

0.0014 

0.400 

1.37 

0.10 

19.10 

-0.26 

0.0005 

0.450 

-0.66 

-0.72 

20.40 

-3.50 

-0.0003 

0.500 

-2.18 

-1.24 

20.73 

-5.60 

-0.0008 

0.600 

-3.86 

-1  .45 

18.61 

-6.67 

-0.0011 

0.700 

-4.09 

-0.82 

13.67 

-4.50 

-0.0004 

0.800 

-3.42 

0.22 

7.54 

-0.65 

0.0012 

0.900 

-2.40 

1.23 

1.83 

3.16 

0.0031 

1  .000 

-1  .53 

1  .85 

-2.29 

5.57 

0.0047 

1  .200 

-1  .20 

1  .50 

-4.29 

4.65 

0.0054 

1.350 

-1.68 

0.35 

-2.12 

0.91 

0.0034 

1  .500 

-2.53 

-0.63 

-0.13 

-2.37 

0.0007 

1  .650 

-2.61 

-0.92 

-0.09 

-3.39 

-0.0007 

1  .800 

-2.12 

-0.61 

-1  .74 

-2.44 

0.0001 

2.100 

-1.13 

-0.10 

-4.1! 

-0.89 

0.0037 

2.400 

-1.16 

-0.74 

-2.07 

-3.07 

0.0016 

2.700 

-0.77 

-0.73 

-1.75 

-3.00 

0.0000 

3.000 

-0.  1 1 

-0.31 

-2.67 

-1.41 

0.0018 

3.600 

-0.16 

-0.61 

-0.77 

-1  .83 

-0.0010 

4.200 

-0.  10 

-0.36 

-0.74 

-0.78 

-0.0008 

4.800 

-0.12 

-0.12 

-0.92 

-0.18 

-0.0014 

5.400 

-0.32 

-0.22 

-0.27 

-0.80 

-0.0022 

6.000 

-0.18 

-0.08 

-0.66 

-0.43 

-0.0024 

Table  4 


WAGNER  FUNCTIONS  FOR  INSTANTANEOUS  AILERON  DEFLECTION  AT  M  =  0.8 


0 

Lift 

W.3(o) 

Pitching 
mome  n  t 

w23(  ) 

Rolling 

moment 

W44(a) 

Bending 
mome  n  t 
W53(o) 

Hinge 
mome  n  t 
W33(o) 

0.0 

0.34  3 

0.470 

0.416 

0.398 

1  .908 

0.  1 

0.429 

0.569 

mssm 

0.502 

1  .699 

0.2 

0.495 

0.643 

■sn 

0.578 

1  .505 

0.3 

0.540 

0.694 

0.669 

0.629 

1.349 

0.4 

0.575 

0.733 

0.708 

0.668 

1  .229 

0.6 

0.630 

0.789 

0.766 

0.728 

1  .090 

0.8 

0.677 

0.832 

0.815 

0.776 

1  .043 

1  .0 

0.715 

0.863 

0.855 

0.815 

1  .031 

1  .2 

0.744 

0.884 

0.882 

0.842 

1  .019 

1  .4 

0.763 

0.895 

0.898 

0.858 

1.001 

1  .6 

0.777 

0.902 

0.908 

0.868 

0.984 

1  .8 

0.787 

0.907 

0.918 

0.875 

0.977 

2.0 

0.797 

0.913 

0.928 

0.882 

0.981 

2.5 

0.819 

0.930 

0.951 

0.900 

0.996 

3.0 

0.836 

0.942 

0.968 

0.914 

0.997 

3.5 

0.853 

0.952 

0.980 

0.926 

0.996 

4.0 

0.869 

0.960 

0.989 

0.936 

0.997 

5.0 

0.900 

0.972 

0.999 

0.949 

0.998 

6.0 

0.928 

0.979 

1  .004 

0.958 

0.998 

7.0 

0.948 

0.984 

1  .004 

0.966 

0.999 

8.0 

0.962 

0.988 

1  .003 

0.975 

0.999 

9.0 

0.971 

0.990 

1  .002 

0.982 

0.999 

10.0 

0.978 

0.992 

1  .000 

0.987 

0.999 

12.0 

0.987 

0.995 

1  .000 

0.992 

1  .000 

14.0 

0.991 

0.996 

1  .000 

0.995 

1  .000 

16.0 

0.993 

0.997 

1  .000 

0.996 

1  .000 

18.0 

0.994 

0.998 

1  .000 

0.997 

1  .000 

20.0 

0.995 

0.998 

1  .000 

0.997 

1  .000 

Steady 

-q;3(o) 

-Q23(0) 

— 

-Q53(°) 

-q;3(°) 

s  ta  te 

=  1.788 

=  1.884 

Hi 

=  5.386 

=  0.0182 

Table  5 


REQUIRED  AILERON  ANGLES,  TAILPLANE  LIFT  AND  RESIDUAL  ROOT  BENDING  MOMENT 


TO  NEUTRALIZE  AN  INCLINED  STEP  DOWN-GUST 


=  0.8,  £ 


tailplane  lift  at  x^/c 


c 

Control 

angles ' 

Tailplane 

lift* 

Bending 

moment* 

Leading 

side 

Trailing 

side 

6t 

0.0 

0.00 

0.00 

-0.000 

0.000 

0.2 

0.23 

0.00 

-0.001 

0.001 

0.5 

1.04 

0.1  1 

-0.004 

0.010 

0.8 

2.44 

0.43 

0.035 

1  . 1 

4.15 

0.99 

-0.026 

0.084 

1  .4 

6.02 

1  .95 

-0.048 

0.169 

1.7 

8.04 

3.29 

-0.077 

0.289 

2.0 

9.93 

4.69 

-0. 107 

0.421 

2.5 

12.15 

6.58 

-0.149 

0.615 

3.0 

13.90 

8.31 

-0.186 

0.782 

3.5 

15.41 

9.94 

-0.220 

0.931 

4.0 

16.72 

11.78 

-0.251 

I  .064 

4.5 

17.75 

13.54 

-0.278 

1.177 

5.0 

18.54 

15.05 

-0 . 300 

1.268 

6.0 

19.40 

18.01 

-0.330 

1  .379 

7.0 

19.51 

19.43 

-0.341 

1  .408 

8.0 

19.36 

19.39 

-0.340 

1  .400 

9.0 

19.21 

19.23 

-0.338 

1.385 

10.0 

19.05 

19.05 

-0.334 

1.370 

12.0 

18.84 

18.82 

-0.329 

1.350 

14.0 

18.81 

18.80 

-0.328 

1.348 

16.0 

18.84 

18.84 

-0.329 

1.351 

20.0 

18.88 

18.88 

-0.329 

1.355 

25.0 

18.87 

18.87 

-0.329 

1  .354 

30.0 

18.87 

18.87 

-0.329 

1 .354 

40.0 

18.86 

18.86 

-0.329 

1  .353 

CO 

18.86 

18.86 

-0.329 

1  .353 

Normalized  with  respect  to  gust  downwash  angle 

Normalized  with  respect  to  steady  magnitude 
on  the  wing  corresponding  to  gust. 
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Table  8 


REQUIRED  AILERON  ANGLES,  TAILPLANE  LIFT  AND  RESIDUAL  ROOT  BENDING  MOMENT 


TO  NEUTRALIZE  AN  INCLINED  (I -COSINE)  DOWN-GUST 


tailplane  lift  at  x  /c 


Control 

angles'*- 

Tailplane 

lift* 

Leading 

side 

6C 

Trailing 

side 

6t 

0.00 

0.00 

-0.000 

0.00 

0.00 

-0.000 

0.01 

0.00 

-0.000 

0.03 

0.01 

-0.000 

0.  10 

0.02 

-0.001 

0.22 

0.06 

-0.001 

0.47 

0.15 

-0.004 

0.87 

0.33 

-0.008 

1  .41 

0.60 

-0.014 

2.11 

0.98 

-0.023 

2.96 

1.47 

-0.034 

3.96 

2.10 

-0.047 

6.32 

3.76 

-0.082 

9.01 

5.92 

-0.125 

11.78 

8.45 

-0.172 

14.33 

11.10 

-0.219 

16.43 

13.61 

-0.261 

18.67 

17.26 

-0.316 

1^.13 

18.81 

-0.333 

19.01 

19.00 

-0.333 

18.88 

18.87 

-0.330 

18.83 

18.85 

-0.329 

18.88 

18.87 

-0.329 

18.87 

18.87 

-0.329 

18.87 

18.87 

-0.329 

18.86 

18.86 

-0.329 

Bending 

moment* 


Normalized  with  respect  to  gust  downwash  angle 

Normalized  with  respect  to  steady  magnitude  on 
the  wing  corresponding  to  the  gust. 
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LIST  OF  SYMBOLS 


instantaneous  bending  moment 

steady  bending  moment  in  second  of  equations  (3-23) 

bending  moment  on  leading,  trailing  side  of  wing  in  equations  (2-12)  to  (2-14) 
maximum  bending  moment 
mean  chord  of  wing 

/H  ;  integer  or  subscript  denoting  force  mode 
integer  or  subscript  denoting  upwash  mode 

Kiissner  function  for  mode  i  in  equations  (2-28)  and  (2-29) 

Kussner  function  for  moment  about  x  =  xT  in  equation  (3-17) 
non-dimensional  loading  due  to  step  gust 
non-dimensional  complex  loading  in  expression  (2-6) 
wavelength  of  gust  in  streamwise,  spanwise  direction 
instantaneous  lift  force 

steady  lift  force  in  first  of  equations  (3-23) 

lift  force  at  location  of  canard,  tailplane 

instantaneous  rolling  moment 

£  for  (1 -cosine)  gust  in  equation  (4-4) 

maximum  rolling  moment 

number  of  spanwise  terms 

Mach  number  of  stream 

number  of  chordwise  terms 

integers  (i  =  1,2,...,N)  (see  Ref  7) 

generalized  coordinate  in  mode  j  in  equations  (3-2)  and  (3-3) 


approximation  to  q^lx)  after  nth  iteration 

generalized  force  indicated  in  equations  (3-20) 

steady  generalized  force  in  second  of  equations  (3-20) 

force  coefficient  in  mode  i  in  equation  (3-4) 

generalized  force  due  to  step  gust  in  equation  (2-25) 

complex  generalized  force  in  equation  (2-7),  Q!.(v)  +  ivQlV(v) 

real  part  of 

semi-span  of  wing 
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LIST  OF  SYMBOLS  (continued) 

time  measured  from  entry  of  wing  into  gust 
time  when  wing  root  leading  edge  enters  gust 
speed  of  stream 

upwash  in  equation  (2-1),  (2-5)  or  (4-1) 

upwash  in  mode  j  (=  1  to  5)  in  equations  (2-8)  to  (2-10)  and 
Wagner  function  for  moment  about  x  =  x^,  in  equation  (3-18) 
Wagner  function  in  equation  (3-1) 

streamwise  coordinate  relative  to  pitching  axis  in  Fig  1 
streamwise  centre  of  pressure 

location  of  gust  front  on  centre-line  at  t  =  0  in  equation  ( 
general  pitching  axis 

location  of  canard  lift,  tailplane  lift 
Length  of  (1-cosine)  gust 
location  of  aileron  hinge  in  Fig  1 
spanwise  coordinate  relative  to  centre-line 
upward  coordinate 

displacement  in  mode  i  (=  I  to  5)  in  equations  (2-8),  (2-10) 
local  incidence  of  (1 -cosine)  gust  in  equation  (4-5) 
final  incidence  of  gust 

instantaneous  aileron  angle  in  streamwise  plane  (radians) 
aileron  angles  on  leading,  trailing  side  of  wing 
maximum  aileron  angle 

non-dimensional  maximum  aileron  rate  in  equation  (3-27) 
spanwise  centre  of  pressure  in  terms  of  y/s 
angle  of  sweepback  of  leading  edge  of  wing 
frequency  parameter,  uc/U 
density  of  stream 

distance  of  travel  in  mean  chords,  Ut/c 
a  relative  to  entry  of  root  leading  edge,  U(t  -  t^/c 
Oq  as  wing  enters  gust  in  equation  (2-21) 

Oq  as  wing  becomes  fully  immersed  (see  Fig  I) 


(2-2) 


—  23) 


and  (2-11) 
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LIST  OF  SYMBOLS  (concluded) 

a  non-dimensional  gust  length,  x  /c 

8  8 

o  the  lesser  of  n  and  c 

m  g 

t  non-dimensional  time,  Ut/c 

i|j  angle  of  inclination  of  gust  front  in  Fig  I 

til  circular  frequency 
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Fig  1 


Fig  1  Planform  with  25%  chord  ailerons  in  inclined  step  gust 
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Fig  3  Transient  rolling  moment  due  to  asymmetric  entry  into  a  step  gust 


Fig  4  Spanwjse  centre  of  pressure  due  to  asymmetric  entry  into  a  step  gust  (M  -  0.8) 
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Fig  5 
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Fig  5  Transient  rolling  moment  from  piston  theory  for  asymmetric  entry  into  a  step  gust 


Fig  6 


Fig  6  Comparison  between  piston  and  lifting-surface  theories  for  step  gust  aerodynamics 
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Leading  side 


Normalized  half  wing  bending  moments  due  to  asymmetric  entry  into  a  step  gust 


Tft  8115? 


1 


Fig  9  Control  motions  to  neutralize  bending  and  rolling  moments  due  to  step  gusts 


Fig  10 


Fig  10  Control  angles  and  tailplane  or  canard  lift  to  neutralize  wing  forces  due  to 
symmetric  gust  entry 
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F*g  1 1  Control  angles  with  tailplane  lift  to  neutralize  wing  forces  due  to  asymmetric  gust  entry 


I 


Fig  12  Asymmetric  control  angles  with  various  locations  of  auxiliary  lift  to 
neutralize  gust  forces  ( 0  =  40°) 
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Fig  14  Banding  moments  with  ailerons  and  tailplane  or  canard  to  neutralize  wing 
forces  due  to  gusts  of  various  inclinations 


Fig  17  Influence  of  gust  length  on  various  force  characteristics 
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Fig  20  Influence  of  gust  length  on  maximum  control  rates 
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